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Abstract 

We calculate the vacuum averages of the energy-momentum tensor 
associated with a massless left-handed spinor fields due to magnetic 
fluxes on idealized cosmic string spacetime. In this analysis three 
distinct configurations of magnetic fields are considered: i) a homoge- 
neous field inside the tube, ii) a, magnetic field proportional to 1/r and 
in) a cylindrical shell with 5- function. In these three cases the axis of 
the infinitely long tubes of radius R coincides with the cosmic string. 
In order to proceed with these calculations we explicitly obtain the Eu- 
clidean Feynman propagators associated with these physical systems. 
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As we shall see, these propagators possess two distinct parts. The 
first are the standard ones, i.e., corresponding to the spinor Green 
functions associated with the massless fermionic fields on the ideal- 
ized cosmic string spacetime with a magnetic flux running through 
the line singularity. The second parts are new, they are due to the 
finite thickness of the radius of the tubes. As we shall see these extra 
parts provide relevant contributions to the vacuum averages of the 
energy-momentum tensor. 



1 Introduction 

Different types of topological defects may have been formed during the phase 
transition in the early universe [1] . Depending on the topology of the vacuum 
manifold they are domain walls, strings, monopoles and textures. Cosmic 
strings have gained some interest recently since they arc considered as a 
good candidate to explain some components of anisotropy on the cosmic 
microwaves background [2], gamma ray bursts [3], gravitational waves [4] 
and highest energy cosmic rays [5]. Also they are thought to be important 
for the structure formation in the universe due to their huge energy density 
per unit length [6]. ^ 

A classical field theory model which presents stringlike solutions is the 
Abehan Higgs model [7] . These solutions, also called "vortices" , correspond 
to infinitely long objects. They have a core radius proportional to the inverse 
of the Higgs mass and magnetic flux tube with radius proportional to the 
inverse of the gauge boson mass. 

Coupling the Abelian Higgs model with gravity, Garfinkle [8] showed that, 
as in fiat spacetime, there exist static cylindrically symmetric solutions rep- 
resenting vortices. He also shown that asymptotically the spacetime around 
the cosmic string is a Minkowiski one minus a wedge. Its core has a non-zero 
thickness, and the magnetic field vanishes outside of it. A complete analysis 
about the structure of the Higgs and magnetic fields near the gauge 
cosmic strings can only be provided numerically. 

The vacuum polarization effects due to a magnetic field confined in a tube 

^For cosmic strings formed at grand unified theory, the energy densities is of order than 
IQi^^Kg/m and their radius 10~^^m. 
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of finite radius in Minkowiski spacetime were first analysed by Serebryanyi 
[9]. A few years later, Guimaraes and Linet [10] and Linet [11], calculated 
tfiese effects for a cfiarged massless scalar and fermionic fields, respectively, 
on a idealized cosmic string spacetime. There a magnetic field running 
through the line singularity was considered. As a consequence, the renor- 
malizcd vacuum expectation vahic associated with the energy-momentum 
tensor, {T^y{x)) Ren.) presents contributions coming from the geometry of the 
spacetime and also the magnetic flux. In order to develop these calculations 
the respective Green functions were obtained. More recently Sriramkumar 
[12] has calculated the vacuum fluctuations of current and energy densi- 
ties for a massless charged scalar field around an idealized cosmic string 
carrying a magnetic flux. There the Green function was obtained taking 
into account the presence of the vector potential in the differential operator, 

= — ieA^, which presents the advantage of calculating the two-point 
function without imposing any boundary condition on the field ^. 

The analysis of the vacuum polarization effects on a massless charged 
scalar field by a magnetic flux confined in a cylindrical tube of finite radius 
in a cosmic string background was performed in [13]. There three distinct 
configurations of magnetic field specified below have been considered. In all of 
them the axis of the infinitely long tube of radius R coincides with the cosmic 
string. Calculating the renormalizcd vacuum expectation value (VEV) of 
the square of the field, ($*(a;)$(a;))Ren, and the energy-momentum tensor, 
{T^{x)) Rem it was observed that these quantities present two contributions 
for each model of magnetic flux. The first are the standard ones due to 
the conical geometry of the spacetime and the magnetic flux. The second 
contributions arc corrections due to the finite thickness of the radius of the 
tube. These extra terms provided relevant contributions for points outside 
the tube. Specifically, for the third model this contribution is a long-range 
effect, i.e., it is as relevant as the standard one up to a distance which exceed 
the radius of the observable universe. 

In this paper we return to this problem and extend the above analysis 
to the charged massless fermionic fields. As we shall see, the renormalized 
VEV of the energy-momentum tensors also present two distinct contribu- 

^Absorbing the gauge field on the definition of the matter fields implies an overall phase 
shift in the Green function. However this procedure provides to {T^v{x)) Ren. the same 
result as obtained considering explicitly the presence of the gauge field in the differential 
operator. 
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tions. Moreover due mainly to interaction of the spin degrees of freedom 
of the fermionic fields with the magnetic fluxes, the corrections due to the 
nonzero thickness of the radius of the tube are composed by two terms. 
Unfortunately it is not possible to provide analytical informations about 
the radial behavior of these corrections, and only numerical analysis can do 
that. Also, additional couplings between the spin with the geometry are also 
present in this system. 

The idealized model for a infinitely long straight static cosmic string 
spacetime can be given in cylindrical coordinates by the line element: 

ds'^ = -dt^ + dr"^ + a^r'^de'^ + dz^ , (1) 

where a = 1 — 4yU is a parameter smaller than unity which codifies the 
presence of a conical two-surface (r, 9) . In fact for a typical Grand Unified 
Theory, a = 1 - O(10~6). 

We shall consider the presence of a magnetic field along the z-direction 
assuming that the field has a finite range in the radial coordinate. We are 
particularly interested in the three models: 

i) H{r) = a^jii ^{R ~ homogeneous field inside; 

ii) H(r) = 2Tvt Rr ^(-^ ~ ^^^'^ proportional to 1/r; 

iii) H{r) = — ^^'^^ ^(r — i?), cylindrical shell, 

where R is the radial extent of the tube, and is the total flux. The ratio of 
the flux to the quantum flux 0o, can be expressed by 0/0o — ^ + 1, where 
N is the integer part and < 7 < 1.^ 

Although the structure of the magnetic field produced by a f/(l)— gauge 
cosmic string cannot be presented by any analytical function, its infiuence 
on the vacuum polarization effects of charged matter fields takes place for 
sure. So in this case, the geometric and magnetic interactions provide con- 
tributions. The relevant physical question is how important they are. In this 
paper, together with previous ones, by Guimaraes and Linet, Linet himself 
and us, we try to answer that question. In our analysis we added an extra 
ingredient on this investigation not yet considered. 

As we have already mentioned, considering a finite thickness for the mag- 
netic fiux, additional contributions to the vacuum averages are obtained. 
Because there is no possibility to analysis these effects for a realistic gauge 
cosmic string, we consider idealized configurations to the magnetic field. This 

"^In this paper we are considering % = G = c = 1. 
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allows us to develop an analytical procedure, and as we shall see, even for 
these configurations, some general informations regarding to the infiuence of 
the nonvanishing thickness of the magnetic flux on the vacuum polarization 
can be extracting. 

This paper is organized as follows. In section 2, we explicitly construct 
the spinor Feynman propagators for three different models of magnetic fields. 
Having these Green functions, in section 3 we calculate the formal expres- 
sion to the respective renormalized vacuum expectation values of the energy- 
momentum tensor ((T^^x)) Ren. associated with the left-handed fermionic fields. 
Because non standard contributions to these VEV do not provide concrete 
informations about their radial behavior, in section 4 we present our numer- 
ical analysis about them for specific values for the parameters a, 7 and A^. 
Finally we leave for section 5 our concluding remarks. 



2 Spinor Feynman Propagator 



The Feynman propagator associated with a charged fermionic field, Sf{x, x'), 
obeys the following differential equation: 

{if + e/l - M) Spix, x') = -^S\x, x')h , (2) 



/-9 

where g = det{—g^i,). The covariant derivative operator reads 

)^ = 4)7^"^^. + r,) , (3) 

e^^^ being the vierbein satisfying the condition ^^a)^'(p)''^^"^'^^^ ~ 



the spin connection given in terms of flat spacetime 7 matrices by 



(4) 



and 



(5) 

In the Appendix we show that if a bispinor Df{x, x') satisfles the differ- 
ential equation 



+e^g'"'A^A, + M^l Df{x, x') = 



'-9 



x')l4 , 



(6) 
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with 



(7) 



TZ being the scalar curvature and the generahzed d'Alembertian operator 
given by 

then the spinor Feynman propagator may be written as 

SFix,x')^{i^ + e4 + M)DFix,x') . (8) 

Now after this brief introduction about the calculation of spinor Feynman 
propagator, let us specialize it to the cosmic string spacetime in the presence 
of a magnetic field along the direction. We shall choose the following base 
tetrad: 

/ 1 \ 

cos^ —sm9/ar 
sin^ cosO/ar 
V 1 / 



(9) 



In order to provide an explicit expression to the differential operators 
above, we shall adopt the same representation as in [15] to the 7— matrices 
in Minkowiski spacetime. 

As to the four vector potential we have 



with 



A^ = (0,0,>l(r),0) , 



(10) 



:ii) 



For the two first models, we can represent the radial function a(r) by: 

a{r) = f{r)Q{R-r) + Q{r -R) , (12) 

with 



f{r) 



r'^/R^^ for the model {%) and 
r/R, for the model {ii). 



(13) 



6 



For the third model, 

a{r) = Q{R-r). (14) 
In this spacetime the only non-zero spin connection is 

r,= '-{l-a)E' (15) 

and for the Christofell symbols we have: 

Tl, = -a'r, Tl, = Tl, = l/r. (16) 

Defining by IC{x) the 4x4 matrix differential operator which acts on 
Df{x,x') in (6), for this physical system we obtain: 

- J^r.i^ - + ^.Mr)de + ^A^(r) + , (17) 

where 
and 

A ^ -d^ + d'^ + -dr + -^^dl + dl . (19) 

r a r 

We can see that the differential operator above explicitly exhibits, besides 
the ordinary d'Alembertian operator on the conical spacetime, four different 
types of interaction terms: (i) the usual charge-magnetic field, (ii) the spin- 
magnetic field, [iii) the spin-geometry and (iv) the spin-charge-geometry. 
All of the last three interactions were absent in the analogous differential 
operator used to define the scalar Green function in [13]. 

Moreover, as we can see the operator )C{x) is diagonal in 2 x 2 blocks. This 
means that the two upper components of the Dirac spinor interact with the 
gravitational and magnetic fields in similar way as the two lower components 
and they do not interact among themselves. 

The system that we want to study consists of masslcss charged fcrmionic 
field in the cosmic string spacetime in the presence of an Abelian magnetic 
field along the direction. Let us chose a left-handed field. In this case the 
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Dirac equation and the equation which defines the spinor Feynman propaga- 
tor reduce themselves to a 2 x 2 matrix differential equation: 



(20) 



where 



(a-^-l)\ 1 



2r 



ar 



a^'^ {de - iA{r)) - a^'^d. 



(21) 



with a^^^ = a.f, a^^^ = a.9 and a^^^ = a.z. 

The Feynman two-components propagator obeys the equation 

1 



and can be given by 



S^{x,x')=iPlG\x,x') , 



where now G^{x,x') obeys the 2x2 matrix differential equation: 



(22) 
(23) 



K{x)G^{x,x') 



'-9 



5\x,x')h , 



(24) 



with 



K{x) = -l^--^(l-a)a^de-eH{r)a^+ ^ 



[l-af 



2ie 



2^2 



1 - a)A{r)a^ + -^A{r)dg + zY^A\r 



Q,2^2 



(25) 



Because of the peculiar diagonal form of the above operator, let us take 
for the following expression: 



G^{x,x') 



G+{x,x') 
G_{x,x') 



(26) 



So the differential equation (24) reduces itself to the two independent ones: 



Q,2^2 



±4^(l-a)A(r)-4^A^ 
a r a r 



G±{x,x') 



1 



-.5\x,x') . 



(27) 
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Due to the cylindrical symmetry of this system each component of the 
Euclidean Green function can be expressed by 



1 pOO POO 

G^ix,x') = E / dk / du;e^^(^-^')e-(---')^±(r,r') .(28) 

(27rj'^ J-00 J-00 

Before to specialize on the specific model, let us write down the non- 

homogcncous differential equation obeyed by the unknown function (/,t(r, r'). 
Substituting (28) into (27) and using the standard representation to the 
delta function in the temporal, angular and 2;-coordinates, we arrive at the 
following differential equation for the unknown function g^ir, r'): 



1 d 
r dr 



±n{l- 2neA) + 



(1 - a^) 



^ e(l - a)A{r) 



(5^±H{r) g^{ry)^—5{T-T') , 



(29) 



where j3'^ = k/ 

It is of our main interest to investigate the vacuum polarization effect for 
external points to the magnetic flux. So, we shall consider solutions of (29) 
with both r and r' greater than R. 

Let us define by 5'jf(r, r') the solution of (29) regular at r — 0, and by 
(7,-^(r, r') the solution that vanishes at infinity. These two solutions must 
satisfy continuity condition at r = r', with their first derivative discontinu- 
ous at this point. Integrating out in the region r < R the inner solutions, 
corresponding to r < r', for the first two models we have: 



for r < it! and 

for < r < r', where 



9n 



4) 



I\,,\{(5r) + E^PR)K\,,\{(3r) 



a 



(30) 



(31) 



(32) 
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with 5 = 1^ = + 7. Hf{r), for i = 1, 2, represents the solution associated 
with the two first models: 

^i^W = jM..(.),|A,,,|(^r^V (33) 



and 



with 



Ht{r) = ^M..(.„|A,(,,| (Cr) , (34) 



n l-a 



_ (^(2n±l) 



(37) 



and 



A,(.) = r ± ^ ■ (38) 



a 2a 

For both functions Hf, ( = -^{6^ + jS'^R^a^f'^ . 

In both cases Mcr^x represents the Whittaker function. The constant 
E'^^-^ (PR) is determined by matching g^"^ and its first derivative at r = i?. So 
we get 

W^^""^ Ht{R)Kl^jm - H'^{R)K\.,\m) ■ ^''^ 

In all the expressions, ly and K^, are the modified Bessel functions. 
The outer solution of (29) is given by 

g%{r, r') = D%^K\^^\{(3r) , for r > r'. (40) 

Now, imposing the boundary conditions on g^"^ and g^^ at r = r', we 
get the following result: 

g^{r,r') = [/|,±|(/3r<) + 4)(/3i?)X|,±|(/3r<)] X|.±|(/3r>) . (41) 

In the above equation (r^) is the larger (smaller) value between r and r' . 
Substituting (41) into (28), developing the sum in n and integrating over k 
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and ijj for the first part, tlie Euclidean Green functions acquire tlie following 
expression: 

^ , e^^^^ e^*^^sinh(7±Mo/a) + sinh[(l -7±)V«] 

' STT^arr' sinh Mo cosh (tio/o;) — cos 

n 

(42) 

with 7=*= = (1 — q;)/2 ^ 7 and 

coshuo = ^ . (43 

2rr 

As we can observe, the first term in these Green functions depends only on 
the conicity parameter and the fractional part of 4>/4>q. Moreover, the second 
part contains information about the radius of the magnetic tube through the 
constants E^y 

As to the third model, the solution to 5'^(r, r') is: 

^±<(r,r')=A±/|.^l(/3r) , (44) 

for r < R, 

^±<(r,r') = S±[7|,^|(/3r)+£;±(/3it:)X|,^l(/3r)] , (45) 
ior R < r < r' and 

gi>ir,r')=D^K^,^^iPr), (46) 

for r > r' . In the above equation e± = l/a{n± {l — a)/2). Again, the coeffi- 
cients A^, and can be determined by imposing boundary conditions 
at r = it! and r — r'. However, due to the expression of the magnetic field 
in this model is concentrate as a 5— function at the cylindrical shell, there 
happens a discontinuity condition in the first derivative of g^^ at r — R. 
The rest being the same. So using these facts we obtain: 

9n{r,r') = -^[/|e^|(/5r<) + E±(/3i?)K|,^l(/5r<)]i^|,^l(/3r>) , (47) 
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where now 



with 



and 

FinaUy, substituting the expression found to above into (28), and 
adopting similar procedure as we did in the two previous cases, we obtain: 

1 g«=FA6» ginh(7-Uo/tt) + sinh[(l — 7)^0/0;] 



Sn'^arr' sinh uq cosh(-Uo/ a) — cos 

X 



00 

Y: e^-^'E^{PR)K^,,0r)K^,,0r') , (49) 

n=— 00 

with 7 = (1 - a)/2. 

Here the first term of the Green function depends only on the conicity 
parameter. 

Although all expressions to and g^ present dependence on the radial 
coordinates r and r', their dependence on r' are implicitly contained in the 
coefficients which multiply the functions of r. These dependences appear as 
consequence of the boundary conditions obeyed by them at r = r'. 



3 Computation of (Too)i?en 



The vacuum expectation value of the energy-momentum tensor associated 
with the system under investigation can be obtained in the book by Grib et al 
[14] , combining respectively the classical expressions to this tensor associated 
with the fermionic field in presence of electromagnetic interaction in a curved 
spacetime. It is given by: 

1 



lim tr 

4 x'^x 



a^{D, - L>,0 - a,{D^ - D^,) 5|(a;, x') , (50) 
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where D^j = Vc. — ieA^r, the bar denotes complex conjugate and cr'' = 

In order to take into account the presence of the three magnetic field 

configurations, given previously, wc write the vector potential in the form 
= (0, 0, ^1^, 0), with a(r) being given by (12) and (13), for the first two 

cases and by (14) for the third case. The spinor Green functions arc expressed 

in terms of the bispinor G^{x,x') given by (26) with G±{x,x') given by (42) 

for the two first models and (49) for the third one. 

For simplicity let us calculate {Tqo{x)) only. The other components of 

this tensor can be obtained by using the conservation condition 

V^{Ti:)nen. = , (51) 

and the vanishing of trace ^ 

{TPn.^. ^ . (52) 

Fortunately only the second order time derivative provides a nonzero contri- 
bution to (Too (a;)). In fact all the other terms go to zero in the coincidence 
hmit and/or after taking the trace over the Pauh matrices. Moreover, be- 
cause the bispinor depends on the time variable with t — t', we finally have: 

(Too(a;)) = - \im tr{d^G{x, x')) 

x'—*x 

= lim tr{d'^GE{x,x')) , (53) 

x'—*x 

where we have made a Wick rotation on the above equation. 

However the calculation of the above expression provides a divergent re- 
sult. In order to obtain a finite and well defined expression, we must apply in 
this calculation some renormalization procedure. The method which we shall 
apply here is the point-splitting renormalization procedure. The basic idea 
of this method consists of subtracting from the Green function all the diver- 
gences which appear in the coincidence limit. In [16] Adler et al observed 
that the singular behavior of the Green function in the coincidence limit 
has the same structure as the Hadamard one. Later Wald [17] introduced a 

'*In general the trace of the renormalized VEV of the energy-momentum tensor is equal 
to j^Tra^ [18]. However, because this spacetime is locally flat and there is no magnetic 
field outside the tube, a-2 = 0. 
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modification to this technique in order to provide the correct resuh for the 
trace anomaly. In this way the renormahzed vacuum expectation value of 
the zero-zero component of the energy-momentum tensor can be given by: 

{Too{x))Ren. = \imtr[d^^GE{x,x')-dlGH{x,x')]. (54) 

x'—^x 

Because this spacetime is locally flat, the Hadamard function coincides with 
the Euclidean Green function in a flat spacetime: 

As it was already mentioned, (42) and (49) present two distinct contribu- 
tions; the first ones contain informations about the geometrical structure of 
the spacetime and the fractional part of the magnetic fiux, and the second, 
the corrections due to the nonzero thickness of the radius of tube. Moreover, 
in the calculation of the VEV, only their first contributions are divergent in 
the coincidence limit, the second ones are finite. Finally, explicitly exhibit- 
ing these remarks, we write down the renormahzed VEV of the zero-zero 
component of the energy-momentum tensor by: 

{Too{x))Ren. = (Tooix)) Reg. + {Too{x))c , (56) 

where 

{Too{x))Reg. = liui \d^^G +{x , x') + d^,G .{x , x') - 2d',G h{x , x')] (57) 

and 

{T^o{x))c^ \i^dlGc{x,x') . (58) 

x'-^x 

Here Gc{x, x') represents the corrections due to the second terms in and 
G- for the three models, as shown next. After some intermediate calculations 
we arrive to the following results: 
i) For the two first models, 

(Too(x))^en. = ^;^^^^[(a2-l)(W + 7) + 1207V-2f -1) 

1 r°'^ °° p 
+ 4^ i 

+El,^{vR/r)Kl^_^{v)] , (59) 
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for i = 1 and 2. 

ii) For the third model, 



{Too{x))Ren. 



— — -(q;2 - 1)(17q;2 + 7) 

1 foo °^ r 

+ 4^/0 ? [E^vR/r)Ki^^{v) 



+E-(vR/r)^_^{v)] . 



(60) 



For all the above expressions the coefficients were given in (39) and (48). 

The new terms obtained in (59) and (60) are consequence of the finite 
thickness of the magnetic flux and will be present even in the absence of 
cosmic string. Unfortunately it is not possible to evaluate these corrections 
analytically. The reason is because the dimcnsionlcss variable of integration 
V appears not only in the argument of the modified Bessel function, but 
mainly because it appears in the order of the Whittaker functions through 
the factors E^. 

Before to provide some qualitative information about the second contri- 
butions of (59) and (60), we would like to make a few comments about our 
results: (i) The first contributions to them depend only on the conicity pa- 
rameter a, and the fractional part of (f)/(f)o, denoted by 7, for the two first 
models, (ii) The second contributions, the corrections, vanish in the limit 
R ^ 0. They also depend on the integer part of 0/0O) ^• 

Some qualitative informations about the behaviors of the second contri- 
butions can be provided: although for the three models these contributions 
present an overall dependence on the radial coordinate, there exist an 
additional dependences in their integrands by the coefficients E"^. However 
quantitative informations about these new behaviors can only be provided 
numerically. So we leave for the next section this analysis. 

4 Numerical Analysis of Corrections 

In this section we shall exhibit the quantitative behavior for corrections to 
the renormalized VEV of the zero-zero components of the energy-momentum 
tensor for the system analysed, considering for the three models specific 
values of the parameters a, 7 and A^. These behavior can only be provided 
by numerical analysis. 
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The numerical method adopted by us to develop this analysis was numer- 
ical routine of MAPLE to evaluate integrals. Our first general conclusion is 
about the contributions presented by each term in the corrections. Although 
the integrals in the variable v present contributions with opposite sign due 
to the coefficient and E^ ^ they do not cancel each other. This fact is 
a consequence, mainly, to the preferential spatial direction dictated by the 
magnetic field on the spinor degrees of freedom. 

Moreover, although the total contributions to the renormahzed VEV of 
the zero-zero component of the energy- momentum tensor, (Too)iien., requires 
in principal an infinite sum of terms, in fact this is not really necessary. The 
sums are very well represented by a few terms only. Expressing this quantity 
as 

(2oo)_Ren. = {Tqo) Reg. (1 " ^'i,a,N{x)) , (61) 

where x — Rjr^ in our graphs we exhibit only the function ^.iv for the three 
models, considering x running in the interval [0.01,0.1]. We also compare 
these functions for the three models. Although, we expect that this function 

goes to zero in the limit as a; — 0, or r ^ cxo, our numerical results are not 
hundred percent trustful for x smaller than 10^^. In fact the behavior of the 
factors E'^^^^i^vx) for the two first models and E^{vx) for the third one for 
small value of x can be obtained by developing a series expansion in power 
of their argument. Our results for them are: 



E, 



{vx) 



2\u±\ 



(62) 



with 

and 
with 



n± 



1 

a 

E^ ~ (vx 



(l-a) 



-S 



e± = - 
a 



n ± 



:i-cx) 



(63) 
(64) 
(65) 



Prom these expressions we can infer that the corrections in fact go to zero 
for large values of the radial distance r, i.e., small value of x. However 
the numerical program provides very irregular diagrams to the integrands of 
(59) and (60) for small values of the dimensionless variable v when x becomes 
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smaller than 10~^, so we can not trust on the evaluations of the integrals for 
these cases. 

In Figs. 1(a) and 1(6), we present the radial dependence of the function 
A^qTV, for = 1 and N — 2, respectively, considering, a — 0.99 and 
7 = 0.02. From both figures we can see that the shapes of the curves are 
very similar, however the larger contributions are due to the third models. 
By our numerical analysis these functions are bigger than unity, indicating 
that the corrections are in fact more relevant than the standard results, 
(Too) Reg., up to a distaucc one thousand times or, even more, larger than the 
radius of the tube. In the Figs. 2(a) and 2(6) we exhibit the logarithmic 
behavior for the same functions for a smaller scale of the variable x. In the 
latter is possible to infer a power behavior for the corrections as expected by 
analytical analysis. 

In Fig. 3, we exhibit the functions „ tv for N = 1 and N = 2. From 
this graph we observe that the correction due to iV = 2 is larger that for 
iV = 1. This result is expected, of course. However we know that the 
values adopted here to the magnetic fluxes are only few times larger than the 
quantum flux, so increasing this magnetic flux for sufficiently large value of 
N the corrections becomes more and more relevant than the standard values 
of (Too) Reg.- 

Finally in Fig. 4, we exhibit the behavior of A^^^^tv varying the parameter 
a for N = 1. From this graphs we observe a prominent increase for this factor 
for values of a bigger than 0.95. Moreover we can observe too that this effect 
becomes more evident in the third model of magnetic field. A similar graph 
has been constructed considering N — 2, however no significant changes were 
found. 

5 Concluding Remarks 

In this paper we have explicitly exhibited the spinor Green functions associ- 
ated with a left-handed charged spinor field on a cosmic string spacetime in 
the presence of external magnetic fiuxes for three different configurations, all 
of them confined inside a long tube of radius R. Having these Green func- 
tions, we calculated the formal expressions to the renormalizcd VEV of the 
zero-zero component of the energy-momentum tensor, {Tqq{x)) R^n- In these 
calculations we observe that two independent contributions were obtained. 
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The first contributions are the standard ones due to the conical geometry 
of the spacetime and the fractional part of <p/<po- They coincide with the 
expression found by Linet in [11]. The second contributions are corrections 
due to the finite thickness of the radius of the tube. Unfortunately, it is not 
possible, by analytical analysis, to provide the complete information about 
the dependence of these parts with the radial coordinate, and only by nu- 
merical analysis we can do this. From this analysis we could observe that 
the corrections are more relevant than the standard results for distances up 
to 10^ bigger than the radius of the tube. Our results also suggest that the 
corrections are almost independent of the specific form of magnetic field, and 
that they increase for higher values of the magnitude of the magnetic field. 
Analysing the behavior of the corrections, ^^,a,N as function of a, we also 
observe that they become more relevant when this parameter becomes closer 
than unity. 

So, from all our numerical analysis we conclude that considering a nonva- 
nishing radial extension to the magnetic fiux in the calculation of the vacuum 
polarization, very important contributions take place. In fact these contri- 
butions become more relevant than the standard ones. The main reason for 
this fact is because no matter how big is the magnitude of the magnetic flux, 
the standard contribution depends only on the fractional part of 0/0o- 

Although the models investigated here in fact do not correspond to the 
realistic configuration associated with a ?7(1)— gauge cosmic string, they pro- 
vide some improvement in the calculation of vacuum polarization effects when 
compared with the ideal case; moreover our results give some general infor- 
mations. They suggest that considering the infiuence of the nonzero thickness 
of the radial extent of the magnetic filed, relevant contributions to the in- 
duced energy densities associated with massless spinor fields are obtained, 
and more, these corrections are almost insensitive to the specific form of the 
magnetic fields. 

Allen et all [20] have investigated the vacuum polarization of a massless 
scalar field on the cosmic string spacetime considering generically the effect 
of the internal structure of the string's core on the metric tensor. In this 
analysis the respective Green function also presents two parts, being one of 
them due to the non-zero core radius of the string. A more realistic treatment 
of the vacuum polarization effect associated with a charged field by a gauge 
cosmic string, must to consider the effect of both non-zero thicknesses for the 
magnetic field and the metric tensor. This analysis deserves to be analyzed 



18 



appropriately in near future. 
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6 Appendix. Square of the Dirac Operator 

In this appendix we want to find the square of the Dirac operator ilp — 
— y + ie/i.. The result presents three distinct contributions. 

a) The gauge field independent term ^: 

i^y^Yv. = + \n , (66) 

being TZ the scalar curvature of the manifold and the generalized d'Alembertian 
operator. 

b) The term linear in the gauge field has two contributions: 

-ie{rV^A,Y + A^rY'^u} , (67) 
which after some intermediate steps can be written as: 

ieig^^^D^A, - W^F,, + 2g^^^A,V,} , (68) 

where 

S'^'^ = ^[7^7l • (69) 

The second term in (68) is a Pauli interaction one. It appears as a conse- 
quence of the square of the Dirac operator. F^y is the usual second-order, 
antisymmetric field-strength tensor. 

c) As to the square gauge field term, it is easily obtained being equal to: 

e^g^^A^A, . (70) 

In our calculations above we have considered that the 7 Dirac matrices 
obey the anticommutator relation {7^,7''} = —2g'^'^. 

^This part is very well known. It is given in the book by Birrell and Davies [19] 
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7 Figure Captions 

Figure 1: These graphs show the behavior of A^^a,N for the three specific 
models of magnetic field configurations with: (a) N — 1, a — 0.99 and 
7 = 0.02. In (6) we change to = 2. 

Figure 2: These graphs present the logarithmic behavior of the previous 
one for A^ = 1 in smaller scale of the variable x. 

Figure 3: This graph shows the behavior of A^^^.a^, for the first model in 
the case N — 1 and N — 2, considering a — 0.99 and 7 = 0.02. 

Figure 4: This diagram presents the behavior of the correction Aj^a,N a-s 
function of a. It was considered x = 10~^, A^ = 1 and S = 0.002. 
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